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Q^ , Abstract. The analysis of the submm anisotropies that will be mapped 

by the forthcoming MAP and PLANCK satellites requires careful fore- 
ground subtraction before measuring CMB fluctuations. Among these, 
the foreground due to IR/submm thermal radiation from dusty sources 
was poorly known until recent observational breakthroughs began un- 
qq ■ veiling the properties of these objects. We hereafter briefly review the 

observational evidence for a strong evolution of IR/submm sources with 
respect to the local universe explored by IRAS. We present the basic 
principles of a new modeling effort where consistent spectral energy dis- 
tributions of galaxies are implemented into the paradigm of hierarchical 
clustering with the fashionable semi-analytic approach. This model pro- 
C*~) ' vides us with specific predictions in IR/submm wavebands, that seem 

, to reproduce the current status of the observations and help assessing 

the capabilities of forthcoming instruments to pursue the exploration of 
the deep universe at IR/submm wavelengths. Finally, the ability of the 
PLANCK High Frequency Instrument all-sky survey to produce a catalogue 
of dusty sources at submm wavelengths is briefly described. 
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1. Introduction 



^ ■ The accurate measurement of the fluctuations of the Cosmic Microwave Back- 

ground down to scales of a few arcmins by the forthcoming satellite missions map 
and PLANCK will require a thorough analysis of all the astrophysical sources of 
submm/mm anisotropies, and a careful separation of the various foreground 
components (Bouchet et al. 1996, Gispert & Bouchet 1997, Tegmark & Efs- 
tathiou, 1996, Tegmark, 1998, Hobson et al. 1998a, Bouchet & Gispert 1999). 
Among them, the foreground due to resolved/unresolved IR/submm galaxies 
that are present at all redshifts on the line of sight was poorly known, until 
observations and analyses in the last three years began unveiling the "optically 
dark" (and infrared-bright) side of galaxy evolution at cosmological distances. 
In parallel to this observational breakthrough, a strong theoretical effort has 
opened up the way to a modelling of these galaxies that was able to implement 
the basic astrophysical processes ruling IR/submm emission in a consistent way. 
As a consequence of this, it is now possible to have a more general view on the 
problems of foreground separation, and on the capabilities of CMB missions to 
put new constraints on the number and properties of these sources. 
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Figure 1. A typical luminosity sequence for nearby spirals, LIRGs and 
ULIRGs, in the spirit of fig. 2 in Sanders & Mirabel (1996). The figure is taken 
from Dcvriendt et al. (1999). The IR luminosities range from 2.1 x 1O 7 L0 for 
the faintest spiral of the sample, to 6 x 1O 1O L0 for M82 (fifth curve from the 
bottom), and 4 x 1O 12 L for the brightest ULIRG of the sample. 

In the local universe, we know from IRAS observations that about 30 % of 
the bolometric luminosity of galaxies is radiated in the IR (Soifer & Neugebauer 
1991). Local galaxies can be classified in a luminosity sequence from spirals (e.g 
the Milky Way - the brightest spirals in the IR have a bar) , and mild starbursts 
(e.g. M82), to the "Luminous Infrared Galaxies" (say, with 1O 11 L < Ljr < 
1O 12 L ), and "Ultra-L uminous Infrared Galaxies" (say, with 10 12 L & < Ljr) 
that radiate more than 95 % of their bolometric luminosity in the IR/submm. 

The IR/submm emission of these sources is due to dust that absorbs UV and 
optical light, and thermally reradiates with a broad spectral energy distribution 
ranging from a few fim to a few mm. Most of the heating is due to young stellar 
populations but, in the faintest objects, the average radiation field due to old 
stellar populations can be the main contributor, and, in the brightest objects 
(especially the ULIRGs) , the question of the fraction of the heating that is due 
to a possible Active Galactic Nucleus is still difficult to assess. However, recent 
work based on ISO observations shows that starbursting still dominates in 80 % 
of ULIRGs, whereas the AGNs power only the brightest objects (Genzel et al. 
1998, Lutz et al. 1998). 

Now, IR/submm observations are beginning to unveil what actually hap- 
pened at higher redshift. The detection of the "Cosmic Infrared Background" 
(hereafter CIRB) at a level twice as high as the "Cosmic Optical Background" 
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(hereafter COB) has shown that about 2/3 of the luminosity budget of galaxies 
is emitted in the IR/submm range (Puget et al. 1996). In the same time, the 
first deep surveys at submm wavelengths have discovered the sources that are 
responsible for the CIRB, with a number density much larger than the usual 
predictions based on our knowledge of the local universe (Smail et al. 1997). 
The optical follow-up of these sources is still in progress, but it appears that 
some (most ?) of them should be the moderate- and high-redshift counterparts 
of the local LIRGs and ULIRGs discovered by the IRAS satellite and thoroughly 
studied by the ISO satellite. 

We shall hereafter focus on these dusty sources with thermal radiation, and 
we refer the reader to the paper by Toffolatti and coworkers (this volume), for 
the description of the foreground due to radiosources that emit free-free and 
synchrotron radiations at larger wavelengths. 
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Figure 2. Observer-frame model spectra of a Ljr = 10 12 i© galaxy at 
increasing redshifts (from top to bottom), for a cosmology with Hq = 50 km 
s _1 Mpc -1 and f2o = 1. The reader is invited to note that the apparent flux 
in the submm range is almost insensitive to redshift, because the shift of the 
100 /jm bump counterbalances distance dimming. 

A consistent approach to the early evolution of galaxies is particularly im- 
portant for any attempt at predicting their submm properties. Three basic 
problems have to be kept in mind, that explain why it is so difficult, starting 
from general ideas about galaxy evolution, to get a correct assessment of the 
number density of faint submm sources and of the level of submm fluctuations 
they generate. 
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First, it is difficult to extrapolate the IR/submm properties of galaxies from 
our knowledge of their optical properties. It is well known that there is no cor- 
relation between the optical and IR fluxes - see e.g. Soifer et al. (1987) for 
an analysis of the statistical properties of the "Bright Galaxy Sample". In- 
terestingly, the galaxies with the highest luminosities also emit most of their 
bolometric luminosity in the IR. If young stars are the dominant source of heat- 
ing, it turns out that the strongest starbursts mainly emit in the IR/submm. 
Fig. 1 shows a sequence of spectral energy distributions for local galaxies with 
various IR luminosities, very much in the spirit of fig. 2 of Sanders & Mirabel 
(1996). The sources have been modelled by Devriendt et al. (1999, see section 
5) under the assumption that starbursts are the dominant source of heating. 
Now, it is known that local LIRGs and ULIRGs are interacting systems and 
mergers (e.g. Sanders & Mirabel 1996). It is consequently plausible that their 
number density should increase with redshift, when more fuel was available for 
star formation and more interactions could trigger it. As a matter of fact, the 
Hubble Deep Field (HDF, Williams et al. 1996) has unveiled a large number 
of irregular /peculiar objects undergoing gravitational interactions (Abraham et 
al. 1996). Such a large number of interacting systems is of course predicted by 
the paradigm of hierarchical clustering, but the quantitative modelling of the 
merging rates of galaxies, and of the influence of merging on star formation is 
highly uncertain. 

Second, we might have kept so far the prejudice that high-redshift galaxies 
have little extinction, simply because their heavy-element abundances are low 
(typically 1/100 to 1/10 of solar at z > 2). However, low abundances do not 
necessarily mean low extinction. For instance, if we assume that dust grains 
have a size distribution similar to the one of our Galaxy (n(a)da oc a -3 ' 5 with 
a-min < a < a max ), and are homogeneously distributed in a region with radius 
R, the optical depth varies as r oc a^n-^i whereas the total dust mass varies 
as Md us t oc a^ ax B?. For a given dust mass and size distribution, there is more 
extinction where grains are small, and close to the heating sources. This is 
probably the reason why Thuan et al. (1998) observed a significant dust emission 
in the extremely metal-poor galaxy SBS0335-052. In this context, modelling 
chemical evolution and transfer is not an easy task. 

Third, distant galaxies are readily observable at submm wavelengths. Fig. 
2 shows model spectra of an ULIRG as it would be observed if placed at different 
redshifts. There is a wavelength range, between ~ 600 |im and ~ 4 mm, in which 
the distance effect is counterbalanced by the "negative k-correction" due to the 
rest-frame emission maximum at ~ 100 /im. In this range, the apparent flux of 
galaxies depends weakly on redshift to the point that, evolution aside, a galaxy 
might be easier to detect at z = 5 than at z = 0.5. The observer-frame submm 
fluxes, faint galaxy counts, and diffuse background of unresolved galaxies are 
consequently very sensitive to the early stages of galaxy evolution. Note that 
this particular wavelength range brackets the maximum of emission of the CMB. 
As a consequence, any uncertainty in the modelling of galaxy evolution at high 
z will strongly reflect on the results of the faint counts of resolved sources, and 
on the fluctuations of the foreground of unresolved sources. 

In sections 2 and 3, we report respectively on the recent observation of 
the CIRB, and the faint submm counts with the ISO satellite and the SCUBA 
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Figure 3. The Cosmic Optical and Infrared Backgrounds. Open triangle, 
solid dots, and open dot : COB obtained from faint counts and compiled by 
Pozzetti et al. (1998). Solid triangle at 0.0912 /im : upper limit by Vogel et 
al. (1995). Solid triangle at 15 : Oliver et al. (1997). Thick solid line 
: CIRB extracted by Fixsen et al. (1998). Thin dotted lines : error bars. 
Broken solid line : CIRB in the cleanest regions of the sky (Puget et al. 1996, 
Guiderdoni et al. 1997). Open squares : upper limits and detections (at 140 
and 240 /un) from Hauser et al. (1998). Solid squares : reanalysis at 100 
and 240 ^m (Lagache et al. 1999) and 3.5 /im (Dwek & Arendt 1998). The 
no-evolution curve (dotted line) is computed from the local IRAS luminosity 
function extrapolated to z = 8, for Hq = 50 km s _1 Mpc -1 and £Iq = 1. The 
thin solid line shows the best model (the so-called "model E" ) in Guiderdoni 
et al. (1998). Dashes and dots give "model A". 

instrument on the James Clerk Maxwell Telescope (see the review by Mann 
and coworkers in this volume). In section 4, we briefly mention the efforts to 
correct the optical surveys for the effect of extinction, that give a lower limit of 
the number of submm sources from the number of sources detected at optical 
wavelengths. In section 5, we resume various attempts developped so far to 
compute consistent optical/IR spectra, and to model IR/submm counts. In 
section 6, we summarize recent developments of the semi-analytic modelling of 
galaxy formation and evolution where the computation of dust extinction and 
emission is explicitly implemented. Finally, in section 7, we sketch an overview 
of the sensitivities of forthcoming instruments that should greatly improve our 
knowledge of IR/submm sources, and we emphasize the capability of the PLANCK 
High Frequency Instrument to get an all-sky survey of bright, dusty sources at 
submm wavelengths. 
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Figure 4. One of the Southern fields observed by ISOPHOT at 175 /jm for 
the FIRBACK program (from Puget et al. 1999). The pixel size is 1.5 arcmin. 
The large-scale fluctuations are due to the presence of Galactic cirrus, and 
24 sources have been identified at 6*175 > 100 mJy. 

2. The Cosmic Infrared Background 

The epoch of galaxy formation can be observed by its imprint on the background 
radiation that is produced by the accumulation of the light of extragalactic 
sources along the line of sight. The direct search for the COB currently gives only 
upper limits. However, estimates of lower limits can be obtained by summing up 
the contributions of faint galaxies. The shallowing of the faint counts obtained 
in the HDF (Williams et al. 1996) suggests that these lower limits are close to 
convergence. 

In the submm range, Puget et al. (1996) have discovered an isotropic com- 
ponent in the FIRAS residuals between 200 fim and 2 mm. This measure was 
confirmed by subsequent work in the cleanest regions of the sky (Guiderdoni et 
al. 1997), and by an independent determination (Fixsen et al. 1998), giving 
a mean value of the background vl v = 1.3 x 10" 5 (Ai o) _0 - 64 ^-B l/ (T d = 18. 5K) 
where A100 is the wavelength in units of 100 /iin. The analysis of the DIRBE 
dark sky has also led to the detection of the isotropic background at 240 and 
140 /iin, and to upper limits at shorter wavelengths down to 2 /im (Schlegel et 
al. 1998, Hauser et al. 1998). Recently, a measure at 3.5 /im was proposed by 
Dwek & Arendt (1998). The results of these analyses seem in good agreement, 
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though the exact level of the background around 140 and 240 fim is still a matter 
of debate. The controversy concerns the correction for the amount of Galactic 
dust in the ionized gas uncorrelated with the HI gas. A new assessment of the 
issue by Lagache et al. (1999) leads to values of the CIRB that are in good 
agreement with the fit of FIR AS data by Fixsen et al. (1998), and to values at 
140 and 240 /im that are lower than in Hauser et al. (1998). Figure 3 displays 
the various determinations. 
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Figure 5. The evolution of the cosmic Star Formation Rate comoving den- 
sity psfr with redshift z. Solid pentagon : local value (Gallego et al. 1995). 
Solid dots and error bars (dotted lines) : CFRS uncorrected for extinction 
(Lilly et al. 1996). Solid dots and error bars (solid lines) : multi-wavelength 
analysis including IR, submm, and radio data (Flores et al. 1999). Open 
squares and error bars (dotted lines) : HDF uncorrected for extinction (Madau 
et al. 1996, 1998). Open squares and error bars (solid lines) : values corrected 
according to Pettini et al. (1998). The corrections derived by Meurer et al. 
(1997) would shift the corrected points upwards by ~0.5 dex. Open triangles 
: (Connolly et al. 1997). The rest-frame UV fluxes are converted into SFRs 
according to Guiderdoni et al. (1998), and a Salpeter IMF. Solid line : best 
model in the latter paper (the so-called "model E"). 

It appears very likely that this isotropic background is the long-sought 
CIRB (Puget et al. 1996, Dwek et al. 1998). As shown in fig. 3, its level is 
about 5-10 times the no-evolution prediction based on the local IR luminosity 
function determined by IRAS. There is about twice as much flux in the CIRB 
than in the COB. If the dust that emits at IR/submm wavelengths is mainly 
heated by young stellar populations, the sum of the fluxes of the CIRB and COB 
gives the level of the Cosmic Background associated with stellar nucleosynthesis 
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(Partridge & Peebles 1967). The bolometric intensity (in W m 2 sr r ) is : 
T f tbol dl eg pzjz = 0) 

b ° l J 47T (1 + ZY 4vT (1 + Z eff ) { ' 

where eboi{t) = r](l + z) 3 pz(t) is the physical emissivity due to young stars 
at cosmic time t, and z e ff is the effective redshift for stellar He and metal 
nucleosynthesis. An approximate census of the local density of heavy elements 
pz(z = 0) ~ 1 x 10 7 Mq Mpc -3 , taking into account the metals in the hot 
gas of rich galaxy clusters (Mushotzky <fe Loewenstein 1997) gives an expected 
bolometric intensity of the background I^i ~ 50(l + z e jj) _1 nW mT 2 sv~ l . This 
value is roughly consistent with the observations for z e ff ~ 1 - 2. 
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Figure 6. Model spectra in the IR and submm, for IR luminosities 10 6 , 
10 , 10 10 , 10 12 , and 10 14 L©. Emissivity index of big grains : m — 2 (dashes), 
m = 1.5 (solid lines), m = 1 (dotted lines). The sequence is drawn from 
Guiderdoni et al. (1998). 

Of course, it is not clear yet whether star formation is responsible for the 
bulk of dust heating, or there is a significant contribution of AGNs. In order to 
address this issue, one has first to identify the sources that are responsible for 
the CIRB. The question of the origin of dust heating in heavily-extinguished 
objects is a difficult one, because both starburst and AGN rejuvenation can be 
fueled by gas inflows triggered by interaction, and IR/submm spectra can be 
very similar if extinction is large. However, according to Genzel et al. (1998), 
the starburst generally contributes to 50-90 % of the heating in local ULIRGs. 
About 80 % of the ULIRGs in the larger local sample of Lutz et al. (1998) are 
dominated by the starburst, but the trend decreases with increasing luminosity, 
and the brightest objects are AGN-dominated. Now, it is very likely that the 
high-redshift counterparts of the local LIRGs and ULIRGs are responsible for 
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the CIRB. However the redshift evolution of the fraction and power of AGNs 
that are harbored in these distant objects is still unknown. 

3. Far infrared galaxies at high redshift 

Various submm surveys have been achieved or are in progress. The FIRBACK 
program is a deep survey of 4 deg 2 at 175 /jm with the ISOPHOT instrument 
aboard ISO. The analysis of about 1/4 of the Southern fields (that is, of 0.25 deg 2 , 
see fig. 4) unveils 24 sources (with a 5a flux limit S v > 100 mJy), corresponding 
to a surface density five times larger than the no-evolution predictions based 
on the local IR luminosity function (Puget et al. 1999). It is likely that we 
are actualy seeing the maximum emission bump at 50-100 [im redshifted at 
cosmological distances, rather than a local population of sources with a very 
cold dust component, which seems to be absent from the shallow ISOPHOT 
survey at 175 //m (Stickel et al. 1998). The total catalogue of the 4 deg 2 will 
include about 275 sources (Dole et al. 1999). The radio and optical follow-up 
for identification is still in progress. This strong evolution is confirmed by the 
other 175 fim deep survey by Kawara et al. (1998). 

The ISOCAM deep surveys at 15 /im also conclude to a significant evolu- 
tion of the sources (Oliver et al. 1997, Aussel et al. 1998, Elbaz et al. 1998, 
1999). Most of the sources identified so far by the optical follow-up have typical 
redshifts z ~ 0.7, and optical colours similar to those of field galaxies, with mor- 
phologies that frequently have signs of interaction. The surveys seem to show 
a population of bright peculiar galaxies, starbursts, LIRGs, and AGNs. The 
observer-frame 15 /im waveband corresponds to rest-frame wavelengths that 
probe the properties of PAH and very small grains, at the depth of the survey. 
The extent to which the 15 /im flux is related to the bulk of IR/submm emission 
produced by star formation is under study. 

Various deep surveys at 850 fim have been achieved with the SCUBA in- 
strument at the JCMT (Smail et al. 1997, Hughes et al. 1998, Barger et al. 
1998, Eales et al. 1998). They also unveil a surface density of sources (with 
S v > 2 mJy) much larger than the no-evolution predictions (by two or three 
orders of magnitude !). The total number of sources discovered in SCUBA deep 
surveys now reaches about 35 (see e.g. Blain et al. 1998) and should rapidly 
increase. The tentative optical identifications seem to show that these objects 
look like distant LIRGs and ULIRGs (Smail et al. 1998, Lilly et al. 1999). In 
the HDF, 4 of the brightest 5 sources seem to lie between redshifts 2 and 4 
(Hughes et al. 1998), but the optical identifications are still a matter of debate 
(Richards, 1998). The source SMM 02399-0136 at z = 2.803, which is gravita- 
tionally amplified by the foreground cluster A370, is clearly an AGN/starburst 
galaxy (Ivison et al. 1998, Prayer et al. 1998). 

4. The optical view and the issue of extinction 

Recent observational breakthroughs have made possible the measurement of the 
Star Formation Rate (SFR) history of the universe from rest-frame UV fluxes of 
moderate- and high-redshift galaxies (Lilly et al. 1996, Madau et al. 1996, 1998, 
Steidel & Hamilton 1993, Steidel et al. 1996, 1999). Since the early versions 
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of the reconstruction of the cosmic SFR density, much work has been done to 
address dust issues. However, a complete assessment of the effect of extinction on 
UV fluxes emitted by young stellar populations, and of the luminosity budget 
of star-forming galaxies is still to come. Dust seems to be present even at 
large redshifts, since the optical spectrum of a gravitationally-lensed galaxy at 
z = 4.92 (Franx et al. 1997) already shows a reddening factor amounting to 
0.1 < E(B - V) < 0.3 (Soifer et al. 1998). 




Figure 7. Spectral energy distribution of an ULIRG, computed by De- 
vriendt et al. (1999). The thin solid line shows the typical spectrum of a 
starburst at age 0.10 Gyr with the same bolometric luminosity. The thick 
solid line and dashes give the spectrum with two different geometries for the 
dust and star distribution : an oblate spheroid (preferred geomety), and a 
screen model. Whereas there is almost no difference in the IR/submm part, 
the predicted optical properties are very sensitive to the dust distribution. 

The cosmic SFR density determined only from the UV fluxes of the Canada- 
France Redshift Survey has been recently revisited with optical, IR, and radio 
observations. The result is an upward correction of the previous values by an 
average factor 2.9 (Flores et al. 1999). At higher redshift, various authors have 
attempted to estimate the extinction correction and to recover the fraction of UV 
starlight absorbed by dust (e.g. Meurer et al. 1997, Pettini et al. 1998). It turns 
out that the observed slope a of the UV spectral energy distribution F\(X) oc \ a 
(say, around 2200 A) is flatter than the standard value Qo — —2.5 computed from 
models of spectrophotometric evolution. The derived extinction corrections are 
large and differ according to the method. For instance, Pettini et al. (1998) and 
coworkers fit a typical extinction curve (the Small Magellanic Cloud one) to the 
observed colors, whereas Meurer et al. (1997) and coworkers use an empirical 
relation between a and the FIR to 2200 A luminosity ratio in local starbursts. 
The former authors derive < E(B — V) >~ 0.09 resulting in a factor 2.7 ab- 
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sorption at 1600 A, whereas the latter derive < E(B — V) >~ 0.30 resulting in 
a factor 10 absorption. This discrepancy suggests sort of a bimodal distribution 
of the young stellar populations : the first method would take into account the 
stars detected in the UV with relatively moderate reddening/extinction, while 
the second one would phenomenologically add the contributions of these "ap- 
parent" stars and of heavily-extinguished stars. Fig. 5 shows the cosmic SFR 
comoving density in the early version (no extinction), and after the work by 
Flores et al. (1999) at z < 1 and the extinction correction derived by Pettini et 
al. (1998) at higher redshift. 

The broad maximum observed at z ~ 1.5 to 3 (see fig. 5) seems to be cor- 
related with the decrease of the cold-gas comoving density in damped Lyman-a 
systems between z = 2 and z = (Lanzetta et al. 1995, Storrie-Lombardi et al. 
1996). These results nicely fit in a view where star formation in bursts triggered 
by interaction/merging consumes and enriches the gas content of galaxies as 
time goes on. It is common wisdom that such a qualitative scenario is expected 
within the paradigm of hierarchical growth of structures. The implementation 
of hierarchical galaxy formation in semi-analytic models confirms this expecta- 
tion (e.g. Baugh et al. 1998, and references therein). The question is to know 
whether the observations of the optically dark side of galaxies could modify this 
view significantly. 
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Figure 8. Spectral energy distribution of a spiral galaxy. The solid line 
shows the fit of the data with a theoretical spectrum computed by Devricndt 
et al. (1999). 
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5. Modelling dust spectra and IR/submm counts 

Various models have been proposed to account for the IR/submm emission of 
galaxies and to predict forthcoming observations. The level of sophistication 
(and complexity) increases from pure luminosity and/or density evolution ex- 
trapolated from the IRAS local luminosity function with (1 + z) n laws, and mod- 
ified black-body spectra, to physically-motivated spectral evolution. 



Name: IRAS 13112+0305 

Type: ULIRG 

L IR = 3.4 10'= L e 




1 10 1 00 1 000 1 4 1 0« 10" 



A (M m ) 

Figure 9. Spectral energy distribution of an ULIRG. The solid line shows 
the fit of the data with a theoretical spectrum computed by Devriendt et al. 
(1999). 

The evolution of the IR/submm luminosities can be computed from the 
usual modelling of spectrophotometric evolution, by implementing the involved 
physical processes (stellar evolutionary tracks and stellar spectra, chemical evo- 
lution, dust formation, dust heating and transfer, dust thermal emission). Dwek 
(1998) tried to explicitly model the processes of dust formation and destruction 
(see references therein for a review of this complicated issue). Most models pre- 
fer to assume simple relations between the dust content and the heavy-element 
abundance of the gas. The simplest assumption is a dust-to-gas ratio that is 
proportional to the heavy-element abundances. 

Guiderdoni et al. (1996, 1997, 1998) proposed a consistent modelling of 
IR/submm spectra that was designed to be subsequently implemented in semi- 
analytic models of galaxy formation and evolution (see section 6). The values of 
the free parameters that appear in this modelling (gas mass and metallicity, ra- 
dius of the gaseous disk) are readily computable in semi-analytic models for the 
overall population of galaxies. The IR/submm spectra of galaxies are computed 
according to Guiderdoni & Rocca-Volmerange (1987), as follows : 

1. follow chemical evolution of the gas ; 
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Figure 10. Predicted counts at 15 /im, 60 /im, 175 fim, and 850 (im for 
models A (dashes and dots) and E (solid line). Model E is normalized to 
reproduce the CIRB. Dotted line : no evolution. Dashes : FIR/submm ther- 
mal sources in Toffolatti et al. (1998). Solid squares : ISO-HDF at 15 /xm 
(Oliver et al. 1997) ; solid hexagon : Rush et al. (1993). Open stars : Faint 
Source Survey (Lonsdale et al. 1990). Open hexagon : ISOPHOT at 175 /zm 
(Kawara et al. 1998) ; triangles : ISOPHOT at 175 ^m (Puget et al. 1999, 
Dole et al. 1999). Open squares : deep SCUBA survey (Smail et al. 1997) ; 
solid dot : SCUBA-HDF (Hughes et al. 1998). 



2. implement extinction curves which depend on metallicity according to ob- 
servations in the Milky Way, the LMC and SMC ; 

3. compute t\ oc Zg as N gas (A\/Av) where Z gas and N gas are the gas metal- 
licity and column density, s = 1.6 for A > 2000 A (and 1.35 below), and 
Ax/Ay is the Milky Way extinction curve ; 

4. assume the so-called "slab" or oblate spheroid geometries where the star 
and dust components are homogeneously mixed with equal height scales. 
The choice of these simple geometries for transfer is motivated by studies 
of nearby samples (Andreani & Franceschini 1996) ; 

5. compute a spectral energy distribution by assuming a mix of various dust 
components (PAH, very small grains, big grains) according to Desert et al. 
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(1990). The contributions are fixed in order to reproduce the observational 
correlation of IRAS colours with total IR luminosity (Soifer & Neugebauer 
1991). 

Fig. 6 gives the luminosity sequence from Guiderdoni et al. (1998). Frances- 
chini et al. (1991, 1994, and other papers of this series) follow the same scheme 
for items 1, 2, 3, 4, and use slightly different IR/submm templates, whereas 
Fall et al. (1996) basically use constant dust-to-metal ratios and black-body 
spectra. Recently, Silva et al. (1998) proposed a more sophisticated treatment 
in which transfer is computed in molecular clouds. The method of Guiderdoni 
et al. (1998) has been extended to obtain far-UV to radio spectra, and to study 
local templates by Devriendt et al. (1999). Fig. 7 shows the predicted opti- 
cal/IR/submm spectrum of an ULIRG, and fig. 8 and 9 display examples of fits 
for observed objects (a spiral galaxy and an ULIRG) from this latter paper. 

These spectra can subsequently be used to model IR/submm counts. The 
simplest idea is to implement luminosity and/or number evolution which are 
parameterized as power laws of (1 + z) (e.g. Blain and Longair 1993a, Pearson 
& Rowan-Robinson 1996 ; see also references in Lonsdale, 1996). These power 
laws are generally derived from fits of the slope of IRAS faint counts (which do 
not probe deeper than z ~ 0.2). Then they are extrapolated up to redshifts of a 
few units. Unfortunately, various analyses of IRAS deep counts yield discrepant 
results at Sqq < 300 mJy, and the amount of evolution is a matter of debate 
(see e.g. Bertin et al. 1997, for a new analysis and discussion). This uncertainty 
increases in the extrapolation at higher z. 

6. Semi analytic modelling 

These classes of models assume that all galaxies form at the same redshift Zf or . 
But the paradigm of the hierarchical growth of structures implies that there is 
no clear-cut redshift Zf or . In this paradigm, the modelling of dissipative and 
non-dissipative processes ruling galaxy formation (halo collapse, cooling, star 
formation, stellar evolution and stellar feedback to the interstellar medium) has 
been achieved at various levels of complexity, in the so-called semi-analytic 
approach which has been successfully applied to the prediction of the statistical 
properties of galaxies (White & Frenk 1991; Lacey & Silk 1991, Kauffmann et 
al. 1993, 1994; Cole et al. 1994; Somerville & Primack 1999 ; see other papers 
of these series). In spite of differences in the details, the conclusions of these 
models in the UV, visible and (stellar) NIR are remarkably similar. 

A first attempt to compute the IR evolution of galaxies with the Press- 
Schechter formalism has been proposed by Blain & Longair (1993a,b), but with 
a crude treatment of dissipative processes. In Guiderdoni et al. (1996, 1997, 
1998), we extend this approach by implementing spectral energy distributions 
in the IR/submm range. As a reference, we take the standard CDM case with 
Hq=50 kms^ 1 Mpc" 1 , = 1, A = and a% = 0.67. We assume a Star 
Formation Rate SFR(t) = M gas /t*, with t* = (3td yn . The efficiency parameter 
1/(3 = 0.01 gives a nice fit of local spirals. The robust result of this type of 
modelling is a cosmic SFR history that is too flat with respect to the data. 

As a phenomenological way of reproducing the steep rise of the cosmic SFR 
history from z = to z = 1, we introduce a "burst" mode of star formation 
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Figure 11. Redshift distributions at 15 /im, 60 fim, 175 fim, and 850 fim 
for models A (dashes and dots) and E (solid line). Model E is normalized to 
reproduce the CIRB. The average values of log 5,, (in Jy) for A log S v = 1 
flux bins are plotted for each curve. For instance, —0.5 stands for 100 mJy 
< S v < 1 Jy. 



involving a mass fraction that increases with z as (1 + z) 4 (similar to the pair 
rate), with ten times higher efficiencies 1/(3 = 0.1. The predicted cosmic SFR 
density is given in fig. 5. The resulting model, called "model A" does not 
predict enough flux at IR/submm wavelengths to reproduce the level of the 
CIRB. This is sort of a minimal model that fits the COB and extrapolates the 
IR/submm fluxes from the optical (see section 4 on extinction). In order to 
increase the IR/submm flux, we have to assume that a small fraction of the gas 
mass (typically less than 10 %) is involved in star formation with a top-heavy 
IMF in heavily-extinguished objects (ULIRG-type galaxies). Then a sequence 
of models is derived. The amount of ULIRG can be normalized e.g on the local 
IRAS luminosity function (to reproduce the IRAS bright counts) and on the level 
of the CIRB. The so-called "model E" normalized to the flux level determined 
in Guiderdoni et al. (1997) (see fig. 3) is hereafter used to predict faint counts 
and redshift distributions. An extension to NIR and visible counts is given in 
Devriendt & Guiderdoni (1999). 

Fig. 10 and 11 give the predicted number counts and redshift distributions 
at 15, 60, 175, and 850 /jm for this model, that have been produced in Guiderdoni 
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L IR =10 12 L bol0 galaxy at z-0.1, 0.5, 1, 3, 5 and 10 
10 N 
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Figure 12. Observer-frame model spectrum of a Ljr — 1O 12 L ULIRG 
at increasing redshifts (from top to bottom), computed by Devriendt et al. 
(1999). The instrumental sensitivities of various past, current, and forthcom- 
ing instruments (ground-based and satellite-borne telescopes) are plotted. 
This figure taken from Guiderdoni et al. (1999) is available upon request to 
guider@iap.fr. Upgrading of the instrumental sensitivities can be achieved 
when it is necessary. 

et al. (1998) before the publication of the ISOPHOT and SCUBA faint counts 
(their fig. 17). The agreement of the predicted number counts with the data 
seems good enough to suggest that these counts do probe the evolving population 
contributing to the CIRB. The model shows that 15 % and 60 % of the CIRB 
respectively at 175 fim and 850 /xm are built up by objects brighter than the 
current limits of ISOPHOT and SCUBA deep fields. The predicted median 
redshift of the ISO-HDF is z ~ 0.8. It increases to z ~ 1.2 for the deep ISOPHOT 
surveys, and to z > 2 for SCUBA, though the latter value seems to be very 
sensitive to the details of evolution. The model by Toffolatti et al. (1998) is also 
shown in fig. 10. It gives more bright sources and less faint sources at submm 
wavelengths. 
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Figure 13. Predicted counts at planck HFI frequencies for models A 
(dashes and dots) and E (solid line). Model E is normalized to reproduce the 
CIRB. Dotted line : no evolution. Dashes : FIR/submm thermal sources in 
Toffolatti et al. (1998). 



7. Future instruments 

Fig. 12 gives the far-UV to submm spectral energy distribution that is typical of 
a Lir = 10 12 Lq ULIRG at various redshifts. This model spectrum is taken from 
the computation of Devriendt et al. (1999). The instrumental sensitivities of 
various past and on-going satellite and ground-based instruments are plotted on 
this diagram : the IRAS Very Faint Source Survey at 60 fim, ISO with ISOCAM sX 
15 fim, and ISOPHOT at 175 /im, the IRAM interferometer at 1.3 mm, SCUBA 
at 450 and 850 jum, and various surveys with the VLA. Forthcoming missions and 
facilities include wire, sirtf, SOFIA, the planck High Frequency Instrument, 
the first Spectral and Photometric Imaging REceiver, and the imaging modes 
of the SUBARU IRCS and VLT VIRMOS instruments. Finally, the capabilities 
of the ngst, LSA/MMA, and Infrared Space Interferometer (darwin) are also 
plotted. 

The final sensitivity of the next-generation instruments observing at IR 
and submm wavelengths (wire, sirtf, SOFIA, planck, first) is going to be 
confusion limited. However, the observation of a large sample of ULIRG-like 
objects in the redshift range 1-5 should be possible. More specifically, the all-sky 
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shallow survey of planck HFI, and the medium-deep survey of first SPIRE 
(to be launched by ESA in 2007), will respectively produce bright (S^ > a few 
100 mJy) and faint (S u > a few 10 mJy) counts that will be complementary. 




z z 



Figure 14. Rcdshift distributions at PLANCK HFI frequencies, for models 
A (dashes and dots) and E (solid line). Model E is normalized to reproduce 
the CIRB. The average values (in Jy) for AlogSV = 1 flux bins are plotted 
for each curve. For instance, —0.5 stands for 100 mJy < S v < 1 Jy. 

Table 1 resumes the various sources of fluctuation in the six bands of 
planck HFI (Bersanelli et al. 1996, and planck HFI Consortium, 1998). 
The confusion limit due to unresolved sources in beam f2 has been roughly 
estimated with the theoretical faint counts from model E, according to the for- 
mula o con f = (/(f Hm S 2 (dN/dS)dS£l) 1/2 . The values a con f and Su m = qa to t 
have been computed iteratively with q = 5. The la total fluctuation is a to t = 
( a ins + a conf + a cir + ^cmb) 1 ^ 2 - However, this does not give a sound estimate of 
what will actually be possible with planck, once proper algorithms of filtering 
and source extraction are implemented. Source extraction can be studied on 
simulated maps. Tegmark & de Oliveira-Costa (1998) showed that (i) atot is 
about 40 to 100 mJy for the HFI frequencies after filtering, (ii) 40 000 and 5000 
sources are expected (at the 5a level in 8 sr) respectively at 857 GHz and 545 
GHz, and (iii) the CMB reconstruction is not jeopardized by the presence of 
point sources at the level predicted by the models. Similar results are obtained 
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by Hobson et al. (1998b) with a maximum-entropy method on mock data (Hob- 
son et al. 1998a). In both estimates, source counts and source confusion are 
based on the predicted counts of Toffolatti et al. (1998), which differ from ours 
: Toffolatti's counts are higher at the bright end, and fainter at the faint end ; 
hence his source confusion is lower. Table 1 gives the number densities expected 
with flux limits of 100 mjy and 500 mjy, according to Guiderdoni et al. (1998). 
The reader is invited to note the strong sensitivity of the counts to the exact 
flux limit. So the expectations of source counts with planck (and first) are 
severely model-dependent. The on-going follow-up of ISOPHOT and SCUBA 
sources will eventually give redshift distributions that will strongly constrain the 
models and help improving the accuracy of the predictions. 

As far as first is concerned, a 10 deg 2 survey with SPIRE will result in 
~ 10 4 sources (first SPIRE Consortium, 1998). The study of the 250/350 
and 350/500 colors are suited to point out sources which are likely to be at 
high redshifts. These sources can be eventually followed at 100 and 170 //m by 
the first Photoconductor Array Camera & Spectrometer and by the FTS mode 
of SPIRE, to get the spectral energy distribution at 200 < A < 600 //m with 
a typical resolution R = A/AA = 20. After a photometric and spectroscopic 
follow-up, the submm observations will readily probe the bulk of (rest-frame 
IR) luminosity associated with star formation. The reconstruction of the cos- 
mic SFR comoving density will thus take into account the correct luminosity 
budget of high-redshift galaxies. However, the spatial resolution of the submm 
instruments will be limited, and only the LSA/MMA should be able to resolve 
the IR/submm sources and study the details of their structure. 



Table 1. Sources of fluctuations in PLANCK HFI wavebands. (1) Cen- 
tral frequency in GHz. (2) Beam full width half maximum in arcmin. 
(3) la instrumental noise for 14 month nominal mission. (4) la fluctu- 
ations due to cirrus at N H i = 1.3 10 20 cm" 2 (level of the cleanest 10 % 
of the sky) . The fluctuations have been estimated following Gautier et 
al. (1992) with P{k) oc k- 2 - 9 and P ,ioo M m = 1-4 lO" 12 ^ 100Mm . (5) la 
CMB fluctuations for AT/T = 10~ 5 . (6) la confusion limit due to IR 
sources in beam fi = 9 2 FWHM . (7) a tot = {^ins+^lonf+^lir+^CMB) 1 ' 2 ■ 
Here a c i r is for Njji = 1.3 10 20 cm -2 . (8) Number counts at a flux 
limit of 100 mJy in 1 sr. (9) Number counts at a flux limit of 500 mjy 
in 1 sr. 
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8. Conclusions 



1. High-redshift galaxies emit much more IR than predictions based on the 
local IR luminosity function, without evolution. The submm counts start 
unveiling the bright end of the population that is responsible for the CIRB. 
The issue of the relative contributions of the starbursts and AGNs to dust 
heating is still unsolved. Local ULIRGs (but the brightest ones) seem to 
be dominated by starburst heating. However the trend at higher redshift 
is unknown. 

2. It is difficult to correct for the influence of dust on the basis of the optical 
spectra alone. Multi-wavelength studies are clearly necessary to address 
the history of the cosmic SFR density. Forthcoming instrument will help 
us greatly improve our knowledge of the optically dark side of galaxy for- 
mation. Next milestones are sirtf, SOFIA, the planck High Frequency 
Instrument, the first Spectral and Photometric Imaging REceiver, and 
the LSA/MMA. 

3. Under the assumption that starburst heating is dominant, simple mod- 
els in the paradigm of hierarchical clustering do reproduce the current 
IR/submm data. These models normalized by means of the current and 
forthcoming counts should help us predict the number of IR/submm sources 
that will be observed by the PLANCK High Frequency Instrument, the con- 
tribution of the unresolved sources to the submm anisotropics, and the 
final strategy for foreground separation and interpretation. Though source 
counts are strongly model-dependent, and only partly constrained by the 
current set of data, the studies so far seem to show that the quality of 
the reconstruction of CMB anisotropies is not severely degraded by the 
presence of foreground point sources. 
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